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Abstract 

We use a new method to estimate with 5% accuracy the contribution of pion and kaon in-flight- 
decays to the dimuon data set acquired with the CDF detector. Based on this improved estimate, 
we show that the total number and the properties of the collected dimuon events are not yet 
accounted for by ordinary sources of dimuons which also include the contributions, as measured in 
the data, of heavy flavor, T, and Drell-Yan production in addition to muons mimicked by hadronic 
punchthrough. 

PACS numbers: 13.85.-t, 14.65.Fy, 13.20.Fc 
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I. INTRODUCTION 



This article presents an improved determination of the composition of a dimuon sample 
recorded in pp collisions at y/s = 1.96 TeV. The data sample consists of events containing 
two central (\rj\ < 0.7) primary (or trigger) muons, each with transverse momentum px > 
3 GeV/c, and with invariant mass larger than 5 GeV/c 2 and smaller than 80 GeV/c 2 . 
The sample may be dominated by real muon pairs due to semileptonic decays of heavy 
flavor, Drell-Yan production and T decays, but also contains events in which one or both 
muons are produced by hadrons that decay in flight or otherwise mimic a muon signal. 
Although the dimuon signature can be a powerful tool with which to search for new physics 
or sources of CP violation, the uncertainty of the in-flight-decay contribution makes the 
precise determination of the fractions of known processes a serious experimental challenge. 
In particular, it remains controversial if muons originating from the decay of objects with 
a lifetime longer than that of heavy-flavored hadrons can be completely accounted for with 
ordinary sources such as in-flight-decays. Earlier and recent studies estimate the fraction 
of this type of event to be negligible |lH3j. Other studies find it significant, suppress it by 
selecting muons produced close to the beamline J], but have estimated its size with a very 
large uncertainty by using Monte Carlo simulations j^J. The present work is basedjm the 



same Monte Carlo simulated samples, and the same analysis methods as Refs. j^, 5], but 
we improve the method to estimate the number of events due to in-flight-decays achieving 
a 5% accuracy. 

Section [III describes the CDF II detector. In Sec. IHIl we review the present experimental 
situation. Sections HV] to |VT] describe the procedure used to tune the simulation and estimate 
the contribution of ordinary sources to events in which muons are produced by objects with 
very long lifetimes. Based on this results, Section IVHI updates the estimate of the rate of 



multi-muon events reported in Ref. J5]. Our conclusions are presented in Sec IVIIIl 

II. CDF II DETECTOR AND TRIGGER 

CDF II is a multipurpose detector, equipped with a charged particle spectrometer and a 
finely segmented calorimeter. In this section, we describe the detector components that are 



relevant to this analysis. The description of these subsystems can be found in Refs. p4l5]. 



8 



Two devices inside the 1.4 T solenoid are used for measuring the momentum of charged 
particles: the silicon vertex detector (SVXII and ISL) and the central tracking chamber 
(COT). The SVXII detector consists of microstrip sensors arranged in six cylindrical shells 
with radii between 1.5 and 10.6 cm, and with a total z coverage 1 of 90 cm. The first SVXII 
layer, also referred to as the LOO detector, is made of single-sided sensors mounted on the 
beryllium beam pipe. The remaining five SVXII layers are made of double-sided sensors 
and are divided into three contiguous five-layer sections along the beam direction z. The 
vertex ^-distribution for pp collisions is approximately described by a Gaussian function 
with a rms of 28 cm. The transverse profile of the Tevatron beam is circular and has a rms 
spread of ~ 25 /im in the horizontal and vertical directions. The SVXII single-hit resolution 
is approximately 11 /im and allows a track impact parameter resolution of approximately 
35 /im, when also including the effect of the beam transverse size. The two additional 
silicon layers of the ISL help to link tracks in the COT to hits in the SVXII. The COT is 
a cylindrical drift chamber containing 96 sense wire layers grouped into eight alternating 
superlayers of axial and stereo wires. Its active volume covers \z\ < 155 cm and 40 to 140 
cm in radius. The transverse momentum resolution of tracks reconstructed using COT hits 
is ct(pt)/pt — 0.0017 [GeV/c] -1 . The trajectory of COT tracks is extrapolated into the 
SVXII detector, and tracks are refitted with additional silicon hits consistent with the track 
extrapolation. 

The central muon detector (CMU) is located around the central electromagnetic and 
hadronic calorimeters, which have a thickness of 5.5 interaction lengths at normal incidence. 
The CMU detector covers a nominal pseudorapidity range \rj\ < 0.63 relative to the center 
of the detector, and is segmented into two barrels of 24 modules, each covering 15° in 0. 
Every module is further segmented into three submodules, each covering 4.2° in <fi and 
consisting of four layers of drift chambers. The smallest drift unit, called a stack, covers 
a 1.2° angle in 0. Adjacent pairs of stacks are combined together into a tower. A track 
segment (hits in two out of four layers of a stack) detected in a tower is referred to as a 

1 In the CDF coordinate system, 9 and <j) are the polar and azimuthal angles of a track, respectively, 
defined with respect to the proton beam direction, z. The pseudorapidity -q is dchned as — In tan(#/2). 
The transverse momentum of a particle is px = p sm(9). The rapidity is defined as y = 1/2 • ln((_E + 
p z )/(E — p z )), where E and p z are the energy and longitudinal momentum of the particle associated with 
the track. 
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CMU stub. A second set of muon drift chambers (CMP) is located behind an additional 
steel absorber of 3.3 interaction lengths. The chambers are 640 cm long and are arranged 
axially to form a box around the central detector. The CMP detector covers a nominal 
pseudorapidity range \rj\ < 0.54 relative to the center of the detector. Muons which produce 
a stub in both the CMU and CMP systems are called CMUP muons. The CMX muon 
detector consists of eight drift chamber layers and scintillation counters positioned behind 
the hadron calorimeter. The CMX detector extends the muon coverage to \rj\ < 1 relative 
to the center of the detector. 

The luminosity is measured using gaseous Cherenkov counters (CLC) that monitor the 
rate of inelastic pp collisions. The inelastic pp cross section at y/s = 1960 GeV is scaled 
from measurements at \fs = 1800 GeV using the calculations in Ref. 16J1. The integrated 

n 

luminosity is determined with a 6% systematic uncertainty [17]. 

CDF uses a three-level trigger system. At Level 1 (LI), data from every beam crossing 
are stored in a pipeline capable of buffering data from 42 beam crossings. The LI trigger 
either rejects events or copies them into one of the four Level 2 (L2) buffers. Events that pass 
the LI and L2 selection criteria are sent to the Level 3 (L3) trigger, a cluster of computers 
running speed-optimized reconstruction code. 

For this study, we select events with two muon candidates identified by the LI and L2 
triggers. The LI trigger uses tracks with > 1.5 GeV/c found by a fast track processor 
(XFT). The XFT examines COT hits from the four axial superlayers and provides r — <fi 
information in azimuthal sections of 1.25°. The XFT passes the track information to a set of 
extrapolation units that determine the CMU towers in which a CMU stub should be found 
if the track is a muon. If a stub is found, a LI CMU primitive is generated. The LI dimuon 
trigger requires at least two CMU primitives, separated by at least two CMU towers. The 
L2 trigger additionally requires that at least one of the muons also has a CMP stub matched 
to an XFT track with pt > 3 GeV/c. All these trigger requirements are emulated by the 
detector simulation on a run-by-run basis. The L3 trigger requires a pair of CMUP muons 
with invariant mass larger than 5 GeV/c 2 , and \5zq\ < 5 cm, where zq is the z coordinate of 
the muon track at its point of closest approach to the beamline in the r — <ft plane. These 
requirements define the dimuon trigger used in this analysis. 
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III. PRESENT UNDERSTANDING OF THE DIMUON SAMPLE COMPOSI- 
TION 



The value of o^-m &->u an d 0"e-s^,c-^, the correlated cross sections for producing pairs of 
central heavy-flavored quarks that decay semileptonically, is derived in Ref. {4] by fitting 



the impact parameter 



18] distribution of the primary muons with the expected shapes from 



all sources believed to be significant: semileptonic heavy flavor decays, prompt quarkonia 
decays, Drell-Yan production, and instrumental backgrounds due to punchthrough of prompt 
or heavy-flavored hadrons which mimic a muon signal [19] . In the following, the sum of these 
processes will be referred to as the prompt plus heavy flavor (P + HF) contribution. The 
notation K puth — > n and n puth — > [i will be used to indicate muon signals mimicked by 
punchthrough of kaons and pions, respectively. In order to properly model the data with 
the templates of the various P + HF sources, the study in Ref. 4| has used strict selection 
criteria, referred to as tight SVX selection in the following, by requiring muon tracks with 
hits in the two innermost layers of the SVX detector, and in at least two of the next four 
outer layers. 

The tight SVX requirements select events in which both muons arise from parent particles 
that have decayed within a distance of ~ 1.5 cm from the pp interaction primary vertex in 
the plane transverse to the beamline. This requirement suppresses the yield of primary 
muons due to in-flight-decays of pions and kaons, in the following referred to as — > fi 
and K l f d — > fi, respectively. This type of contribution to the dimuon dataset prior to any 
SVX requirement was considered negligible in previous and recent [3[ studies by the 

CDF and DO collaborations. 

As shown by Fig. [TJ the tight SVX sample is well modeled by fits using the prompt and 
heavy flavor contributions [4] . The sample composition determined by the fit and corrected 
for the appropriate efficiency of the tight SVX requirements 2 is listed in the first two columns 
of Table E 

The difference between the total number of dimuons and the P+HF component indicates 
the presence of an important source of dimuons produced beyond 1.5 cm which is suppressed 



2 The efficiency of the tight SVX selection has been measured [5j to be 0.257 ± 0.004 for prompt dimuons 
and 0.237 ± O.OOf for dimuons produced by heavy flavor decays by using control samples of data from 



various sources [ J/ip — > fJ. + fJ. , — > K ± , B — > /j,D°, and T — > ). 
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TABLE I: Number of events attributed to the different dimuon sources by the fit to the muon 
impact-parameter distribution. The fit parameters BB, CC, and PP represent the bb, cc, and 
prompt dimuon contributions, respectively. The component BC represents events containing b 
and c quarks. The fit parameter BP (CP) estimates the number of events in which there is only 
one b (c) quark in the detector acceptance and the second muon is produced by prompt hadrons 
in the recoiling jet that mimic a muon signal. Real muons are muons from semileptonic decay of 
heavy flavors, Drell-Yan production or quarkonia decays. The data correspond to an integrated 
luminosity of 742 pb _1 . The dimuon data set consists of 743006 events. 



Component No. of Events 
BB 230308 ± 2861 



CC 



PP 



BP 
CP 
BC 

P + HF 



103198 ± 6603 



161696 ± 2533 



43096 ± 3087 
41582 ± 5427 
9135 ± 2924 
589015 ± 5074 



No. of real fi — fi 
R bb x BB 

Rcc ^ CC 

T = 51680 ± 649 
+ DY = 54200 ± 5420 



No. and type of misidentified fi 
R bb x [7902 KP uth ft 

+8145 irP uth -> fi } 
R cc x [17546 KP uth -»• fi 
+ 9535 irP uth ->• fi] 
4400 KP uth -> fi KP uth -> fi 



+30000 Tt puth 



fl 7P 



puth 



+ 23000 KP uth -»■ fi irP uth 



/' 

> fi 



11909 K puth 
16447 ETf"^ 



fi + 29253 vrf"*' 1 -»■ fi 
fi + 35275 vrf"*' 1 -> /x 



by the tight SVX requirements. Because unnoticed by previous experiments, this source was 
whimsically referred to as the ghost contribution. 

The relative size of the ghost and P + HF contributions depends upon the type of SVX 
requirement applied to the trigger muons. Reference shows that neglecting the presence of 



ghost events affected previous measurements of cr^^ ^^ [I Q and of x [2] at the Tevatron. 
Finally, the ghost sample is shown to be the source of the dimuon invariant mass discrepancy 
observed in Ref. j^ . 

Reference has studied a number of potential sources of muons originating beyond 
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FIG. 1: The projection of the two-dimensional impact parameter distribution of muon pairs onto 
one of the two axes is compared to the fit result (histogram). 



the beam pipe. Contrary to what assumed by previous experiments, the one source found 
to contribute significantly arises from in-flight-decays of pions and kaons. Based upon a 
generic QCD simulation, that study estimates a contribution of 57000 events. A smaller 
contribution (12052 ± 466 events) from K s and hyperon decays in which the punchthrough 
of a hadronic prong mimics a muon signal was estimated using the data. Secondary inelastic 
interactions in the tracking volume were found to be a negligible source of ghost events. 
The final estimate of the size of possible sources of ghost events underpredicts the observed 
number by approximately a factor of two (154000 observed and 69000 accounted for), but 
the difference was not considered significant because of the simulation uncertainty. 

The present study uses events selected with the tight SVX requirements to tune the QCD 
simulation. Since these data are well modeled by the impact parameter templates of the 
P + HF components, misidentified muons can only arise from the punchthrough of prompt 
hadrons or hadrons produced by heavy-flavor decays. The numbers of misidentified muons 
in the data are derived by subtracting the expected number of real muons, listed in the third 
column of Table [H from the corresponding components in the second column. We then com- 
pare these differences to the rate of K puth — > \x and r K vuth — > \x misidentifications predicted 
by the simulation, and listed on the fourth column of the same table. The simulation is 
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tuned by adjusting the predicted rate of pions and kaons to reproduce the observed number 
of muon misidentifications. Then, the tuned simulation is used to predict the number of 
muons due to in-flight-decays with 5% accuracy. 



IV. RATES OF MISIDENTIFIED MUONS IN THE DATA AND SIMULATION 

We make use of three different samples of simulated events generated with the herwig 
parton-shower Monte-Carlo program [22j, the settings of which are described in Appendix A 



of Ref. [|]. We use option 1500 of the herwig program to generate final states produced 
by hard scattering of partons with transverse momentum larger than 3 GeV/c (sample 
A=generic QCD). Hadrons with heavy flavors are subsequently decayed using the evtgen 
Monte Carlo program 23|. The detector response to particles produced by the above gen- 

1 the CDF II detector simulation that in turn is based on the geant 



erators is modeled wit 
Monte Carlo program 



24j . The values of the heavy flavor cross sections predicted by the gen- 



erator are scaled to the measured values cr^ ^ = 1549 ±133 pb and CTc-^s-j^ = 624 ±104 



pb {4]. The next simulated sample (sample B=single b + single c) is extracted from A by 
requiring the presence of at least a trigger muon generated from heavy-flavor semileptonic 
decays. The simulated sample C=bb ± cc is extracted from B by requiring the presence of 
at least two trigger muons generated from heavy flavor decays. This sample has been used 
to construct impact parameter templates and estimate kinematic acceptances in Ref. 4]. 

In the various simulations, we evaluate the number of dimuons from heavy flavor decays 
and the number of pairs of tracks of different type that pass the same kinematic selection. 
The ratios of these numbers are listed in Tables [III to IIV1 The rate of pairs of tracks of 
different type predicted by the simulation are normalized to the data by multiplying these 
ratios by the number of dimuons from bb or cc production observed in the data. 

TABLE II: Ratio of the numbers of tttt, KK, and Kit pairs to that of primary dimuons from bb 
decays (221096 pairs) in the generic QCD simulation (sample A). 



Process 


Rkk 


Rk-k 


R7T7V 


generic QCD 


867 


8935 


22913 
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TABLE III: Ratio of the numbers of \i — K{tt) combinations to that of primary dimuons from 
bb and cc production in the single-6 and single-c simulated samples (221096 and 83590 dimuons, 
respectively) . 

Process Rk R-k 

single b 11.1 54.4 

single c 40.7 173.7 



TABLE IV: Ratio of the numbers of ^ — K{tt) combinations to that of primary dimuons from heavy 
flavor production in the bb and cc simulated samples (221096 and 83590 dimuons, respectively). 

Process Rk Rw 

bb 7.4 15.2 

cc 43.5 46.9 



The probability PfcL) that a kaon (pion) is not contained by the calorimeter and mimics a 
muon signal has been measured in Ref. [4j by using kaons and pions from — > ^D with 
D° — > K + n~ decays. The probability that kaon (pion) in-flight-decays mimic a trigger muon, 
Pk(L\, has been derived in Ref. by using the simulated sample C. These probabilities 
depend on the particle transverse momentum. Table |V] lists the average probabilities that 
kaons (pions) mimic a primary muon when applying the /tf and Pj^u) probabilities to 
simulated kaon (pion) tracks with p T > 3 GeV/c and \r]\ < 0.7. 



TABLE V: Average probabilities (%) that punchthroughs or in-flight decays result into a primary 
muon. The pt distribution of kaons and pions in the different simulations are almost indistinguish- 
able. 



p puth 

<, r K > 


< pr th > 


<P]( d > 


< pi/ d > 


0.483 ± 0.003 


0.243 ± .0.004 


0.345 ± 0.005 


0.0727 ± 0.0016 



By weighting simulated pion (kaon) tracks that pass the muon kinematic selection with 
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the corresponding Pj^* \ probability, we obtain the prediction of misidentified primary muons 
for the various P + HF components that is listed in the fourth column of Table [H The third 
column of the same table lists the number of real muons for the various P+HF contributions. 
The sum of real plus misidentified muon pairs is in general agreement with the data listed 
in the second column of the table. Therefore, it is reasonable to use the observed rate of 
dimuons, the knowledge of the fraction of real dimuons due to semileptonic decay of heavy 
flavors, Drell-Yan or T mesons, and the knowledge of the P^L\ probabilities to normalize 
the absolute yields of pions and kaons predicted by the simulation. The simulation fitted to 
the data is then used to predict the rate of events due to in-flight-decay misidentifications 
by weighting simulated tracks with the probabilities, the average of which is listed in 

Table [V] In addition, the total rate of K — > fi — K puth — > fi + K l ^ d — > fi misidentifications 
predicted by the simulation can be further constrained with data. This is done in the next 
section by using the number of primary muons due to misidentification of K*°, K^, and 
Kg decays. 

We first describe the evaluation of the content of real muons in the various P + HF 
components and the function used to fit the simulation to the data. Reference (4] estimates 
that the fraction R bb = 0.96 ± 0.04 of the BB component is due to real muons from b- 
quark semileptonic decays whereas the remaining 4% is due to muons mimicked by the 
punchthrough of hadrons produced by heavy flavor decays. Similarly, the fraction R cc = 
0.81 ± 0.09 of the CC component is due to real muons from c-quark semileptonic decays 
whereas the remaining 19% is due to muons mimicked by the punchthrough of hadrons 
produced by heavy flavor decays. The uncertainty of the fraction of real muons due to 
bb (cc) production is accounted for by multiplying R bb (R cc ) by the fit parameter f bb (f cc ) 
constrained to 1 with a 4% (11%) Gaussian error. 

The number of T mesons contributing to the PP component (T = 51680 ± 649 candi- 
dates) has been determined in Ref. [4] by fitting the dimuon invariant mass spectrum with 
three Gaussian functions to model the signal and a straight line to model the combinatorial 
background. The Drell-Yan contribution is evaluated as DY = T x (Joy/ax- The cross 
section cr^y in the 5 — 80 GeV/c 2 mass range is evaluated with a NLO calculation 25(, and 
we use the measured value of o~y 26]. The ratio aoy/cr is 1.05 with a 10% error mostly 
due to the measurement in Ref. 26(. To account for the uncertainty, we weight the DY 
contribution with the fit parameter fd y constrained to 1 with a 10% Gaussian error. 
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The magnitude of the BP (CP) component, predicted with the single-5 (single-c) simula- 
tion, with respect to that of the BB (CC) contribution depends on the ratio of NLO to LO 
terms evaluated by the herwig generator. Because of the dependence on the renormaliza- 
tion and factorization scales, the uncertainty of the single-6 (c) cross section to that of the 



bb (cc) cross section is estimated 



27[ to be ~ 20 (30)% 3 . We account for this uncertainty 



by weighting the rate of pion and kaon tracks predicted by the single-6 (single-c) simulation 
with the additional fit parameter f s b (f sc ) constrained to 1 with a 20% (30%) Gaussian error. 

The simulation prediction of the number of muons mimicked by the punchthrough of 
pions (kaons) is weighted with the fit free parameter f n (fx)- These fit parameters provide 
the absolute normalization of the pion (kaon) rate predicted by the simulation including the 
uncertainties of the punchthrough probabilities. 



V. MEASUREMENT OF THE K -> \i CONTRIBUTION 

The small rate of K ->- fi = K puth ->• // + K ifd — > /i misidentifications is measured using 
a higher statistics sample of dimuon events corresponding to an integrated luminosity of 3.9 
fb _1 . The number of K — > fi misidentification, Nk, is derived from N^*o, the number of 
identified K*° — > K + n~ decays with K + — > fi + (and charge-conjugate states). The number 
N K *o is related to N K by 

N K *o = N K -e -R(K*°), 

where R(K*°) is the fraction of kaons that result from K*° — > K + ix~ decays and eo is the 
efficiency to reconstruct the pion. 

We also select K$ — > -k + -k^ with 7T — > jj, candidates and reconstruct K*^ — > decays. 
The number of K*^ is related to that of Kg by 



N K *± = N K o ■ e x ■ R(K 

where R^*^) is the fraction of K$ resulting from K*^ — > decays and €\ is the 

efficiency to reconstruct the additional pion. We use isospin invariance to set R^K*^) = 
R(K*°). Since the additional pion used to search for the and K*° candidates is selected 



However, the study in Ref. 28( shows that the HERWIG generator predicts the observed single and 



correlated heavy-flavor cross sections to better than 10%. 
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with the same kinematic requirements, we set eo = e±. It follows that 



N, 



K 



N K o/N K ,± x N K *o. 



We search for K*° decays by combining primary muons, assumed to be kaons, with all 
opposite charge tracks, assumed to be pions, with pt > 0.5 GeV/c and in an angular cone 
with cos# > 0.6 around the direction of the primary muon. We require tracks with at least 
10 axial and 10 stereo COT hits. We constrain the pair to arise from a common three- 
dimensional point, and reject combinations if the probability of the vertex-constrained fit is 
smaller than 0.001. The invariant mass spectrum of the selected K*° — > K + ir~ candidates 
is shown in Fig. [2j We fit the invariant mass distribution with a Breit-Wigner function 
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FIG. 2: Invariant mass distribution of K*° — > K + ir candidates passing our selection criteria. The 
line represents the fit described in the text. 



smeared by the detector mass resolution to model the signal. We fix the mass and width of 
the Breit-Wigner function to 896 and 51 MeV/c 2 [29], respectively 4 . We use a fourth order 
polynomial to model the combinatorial background under the signal, and the fitted range 
of invariant mass is conveniently chosen to yield a fit with 50% probability. The size of the 
signal is not affected by the arbitrary choice of the function used to model the combinatorial 



4 The mass resolution due to the track reconstruction in simulated events which include kaon in-flight- 
decays is 4.9 MeV/c 2 , and is negligible compared to the resonance width. 
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background or of the fitted mass range, and is solely determined by the accurate knowledge 
of the signal shape. The fit yields N K *o = 87471 ± 2217 K*° mesons. 

We search for Kg — > 7r + 7r~ with a tt — > fi misidentification by combining primary muons 
with tracks passing the same requirements as those used in the K*° search. In this case, 
both tracks are assumed to be pions. As in the previous case, we select pairs consistent with 
arising from a common three-dimensional vertex. We take advantage of the K s long lifetime 
to suppress the combinatorial background. We further require that the distance between the 
K s vertex and the event primary vertex, corrected by the Kg Lorentz boost, corresponds to 
ct > 0.1 cm. The invariant mass of the K s candidates is shown in Fig. [3J 
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FIG. 3: Invariant mass distribution of K s candidates passing our selection criteria. The line 
represents the fit described in the text. 



We fit the signal with two Gaussian functions and the combinatorial background with a 
straight line in the mass range 0.4 — 0.6 GeV/c 2 . Having fixed the peak of the Gaussian 



functions at 0.497 GeV/c 2 29j, the fit returns an averaged a of 8.4 MeV/c 2 , consistent with 
what is expected from simulated events 5 , and a signal of 32445 ± 421 K s mesons in the 
mass range 0.474 - 0.522 GeV/c 2 . 



5 Because of the Kg long decay path, reconstructed track segments may be shorter than the available 
tracking detector length. When Kg mesons decay before entering the COT volume, the mass resolution 
is 4 MeV/c 2 . 
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We search for K by combining K s candidates with mass between 0.474 and 0.522 
GeV/c 2 and ct > 0.1 cm with any additional track, assumed to be a pion, that pass the 
same selection as pion tracks used to find K*° candidates. We constrain the Kg mass to 
0.497 GeV/c 2 and require that the K s candidate and the pion track are consistent with 
arising from a common three-dimensional vertex. The invariant mass distribution of K*^ 
candidates is shown in Fig. HJ 
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FIG. 4: Invariant mass distribution of K*^ — > K s tt^ candidates passing our selection criteria. The 
line represents the fit described in the text. 



We fit the invariant mass distribution with a Breit-Wigner function to model the signal 
and a fourth order polynomial to model the combinatorial background. We fix the mass 
and width of the Breit-Wigner function to 892 and 51 MeV/c 2 |29|, respectively 6 . The fit 
returns a signal of 3326 ± 246 K*^ mesons. 

The signals obtained by analyzing the 3.9 fb _1 sample are rescaled to estimate the number 
of K — > fi misidentifications present in the 742 pb _1 dataset. After rescaling, we obtain 
Nk = N k o/Nk*± x Nk*o = 164769 ± 13067. This number is used to constrain the total 
number iVJ£ m of K — > misidentifications predicted by the simulation. 



6 In simulated events, when constraining the Kg mass to the PDG value, the mass-constrained Kg mo- 
mentum is measured as accurately as that of a track corresponding to a K — > /j, decay. The resulting K* ± 
mass resolution is approximately 5 MeV/c 2 . 
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VI. FIT OF THE SIMULATION PREDICTION TO THE DATA 



We fit the simulation prediction with a ^-minimization method 30] . The \ 2 function is 
defined as 

5 

X 2 = ^2(D[i] - P[i]) 2 /ED[i} 2 + (N s ™ - 164769) 2 /13067 2 , 

i=i 

where D[i] and ED[i] are the size and error of the component listed in the second column 
and i-th row of Table [H The term P[i] is the sum of the real muon contribution and the 
punchthrough contribution predicted by the simulation in the third and forth columns of the 
same table, respectively. These contributions are weighted with the fit parameters described 
in the previous section, and the terms read 



P[l) = f bb (0.96 ■ 230308 + f k 7902 + / w 8145), 

P[2] = f cc (0.81 ■ 103198 + f k 17546 + / w 9535), 

P[3] = T + f dy DY + f K 4400 + f 2 30000 + f K U 23000, 

P[4] = U (fK 11909 + U 29253), and 

P[5] = f, c (f K 16447 + U 35275). 

In addition, the sum Yli=i is constrained to the observed number of P + HF — BC 
events within its error. The fit results are shown in Tables [VI] to IVIIIl In Table IVI[ the fit 
parameters that tune the various cross sections predicted by the herwig generator are very 
close to their nominal values indicating that the default simulation provides a quite accurate 
modeling of the data. 

The fit returns 163501 K — > // candidates (164769 ± 13067 are measured in the data), 
51% of which are due to punchthrough and 49% to in-flight-decays. We verify this result by 
measuring the fraction of K — > /i decays in identified K*° — > K + tt~ decays that pass the tight 
SVX requirements. The efficiencies of the tight SVX requirement applied to primary muons 
are 0.356 ± 0.002 for muons due to punchthrough of prompt and heavy-flavored hadrons 
and 0.166 ± 0.005 for muons arising from in-flight decays 7 . Based on the kaon composition 

7 In Ref. [5], which uses 0.74 fb _1 of data, these efficiencies have been measured to be 0.45 and 0.21, 
respectively. In 3.9 fb _1 of data, by using T candidates, we measure a smaller efficiency of the tight SVX 
selection. The efficiency loss comes from periods of data taking in which the pedestals of the LOO channels 
were miscalibrated. 
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TABLE VI: Parameter values returned by the fit described in the text. The fit yields x 2 = 2.5 for 
5 DOF. 



f bb 0.97 ± 0.01 

fee 0.95 ± 0.04 

fdy 1.01 ± 0.09 

U 0.97 ± 0.08 

f K 1.01 ± 0.08 

/ s6 1.05 ± 0.08 

f sc 0.82 ±0.10 



TABLE VII: Parameter correlation coefficients returned by the fit. 



Fit parameter 


fbb 


fee 


fdy 


U 


fK 


fsb 


fee 


-0.11 












fdy 


0.17 


0.10 










u 


-0.16 


-0.10 


-0.84 








fx 


-0.08 


-0.17 


0.71 


-0.48 






fsb 


-0.07 


-0.03 


0.50 


-0.51 


-0.01 




fsc 


-0.02 


-0.21 


0.36 


-0.35 


-0.12 


0.15 



returned by the fit, we estimate the efficiency of the tight SVX requirement applied to 
K — > jj, misidentifications to be 0.263 ± 0.008, where the error includes the uncertainty of 
the efficiencies and that of the kaon composition returned by the fit. Figure [5] shows the 
invariant mass distribution of K*° candidates after applying the tight SVX requirement. 
We fit the invariant mass distribution with the same function used to fit the K*° mass 
distribution in Fig. |2j The fit returns 22689 ± 985 K*° candidates to be compared with 
87741 ±2219 K*° candidates before applying the tight SVX requirement. The resulting 
efficiency of the tight SVX requirement is 0.253 ± 0.013, in agreement with what expected 
(0.263 ± 0.008) using the composition of the kaon sample returned by the fit. 
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TABLE VIII: Number of events due to different production mechanisms are compared to the result 
of the present fit. 



Component 


No. of Events 


Fit result 


BB 


230308 ± 2861 


230607 


CC 


103198 ± 6603 


104463 


PP 


161696 ± 2533 


161387 


BP 


43096 ± 3087 


42490 


CP 


41582 ± 5427 


41822 


P + HF 


589015 ± 5074 


589905 




0.7 0.8 0.9 1 1.1 1.2 

M (GeV/c 2 ) 

FIG. 5: Invariant mass distribution of K*° candidates in which the K — >■ \i misidentification passes 
the tight SVX requirement. The line represents the fit described in the text. 

Continuing with the analysis of the results returned by the fit, the total number of 
tt — y fi misidentifications is 240915, 64% of which are due to punchthrough and 36% to 
in-flight-decays. The fractional composition of the K — >■ fi and tt — > fi misidentifications is 
summarized in Table IIXI 

The total fraction of misidentified muons in the dataset is 27%. The number of misidenti- 
fied muons due in-flight-decays of pions and kaons (ghost events) is 113613 ±5332. Since the 
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TABLE IX: Contributions (%) of various processes to pion or kaon misidentifications. 



Type 


R puth _^ p 




All sample 


51 


64 


generic QCD 


20 


33 


single b + bb + 






single c + cc 


31 


31 


Type 






All sample 


49 


36 


generic QCD 


27 


27 


single b + bb + 






single c + cc 


22 
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number of muons from in-flight-decays is derived from that of muons mimicked by hadron 
punchthrough using the fake probabilities listed in Table |Vl the uncertainty of these proba- 
bilities yields an additional error of 3845 events. 

After adding the 12052±466 events from Kg and hyperon decays, we predict 125665±5351 
ghost events, whereas the dimuon dataset contains 153991 ± 5074 events of this type. The 
number of unaccounted events (28326 ± 7374) is (12.8 ± 3.2)% of the bb production and 
(18.3 ± 4.7)% of the ghost sample. 

VII. REVISED ESTIMATE OF THE RATE OF ADDITIONAL REAL MUONS IN 
GHOST EVENTS 

As a cross-check of its bb content, Reference Q has investigated the rate of sequential 
semileptonic decays of single b quarks in the dimuon sample. We provide here a summary of 
that study and its conclusions. That study searches for additional muons with px > 2 GeV/c 
and \q\ < 1.1 in a dimuon sample corresponding to an integrated luminosity of 1426 pb _1 . 
The sample of 1426571 events contains 1131090 ±9271 P + HF events in which both muons 
originate inside the beam pipe and 295481 ± 9271 ghost events in which at least one muon 
is produced outside. The study selects pairs of primary and additional muons with opposite 
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charge {OS) and invariant mass smaller than 5 GeV/c 2 . 

In the case of Drell-Yan or quarkonia production, which was not simulated, the rate 
of same-charge pairs (SS) is a measure of the fake muon contribution since misidentified 
muons arise from the underlying event which has no charge correlation with primary muons. 
The rate of additional muons mimicked by hadronic punchthrough is also estimated with 
a probability per track derived by using kaons and pions from D*^ — > 7r ± D° with D° — >■ 
K + n~ decays. This misidentification probability is approximately ten times larger than 
that for primary muons that have to penetrate twice as many interaction lengths 8 . The 
punchthrough probabilities for pions and kaons differ by a factor of two. In addition, in 
simulated events due to heavy flavor production, the pion to kaon ratio depends on the 
invariant mass and the charge of the muon-hadron pairs. Therefore, for P + HF events, 
the rate of OS — SS pairs is compared to that predicted by the heavy flavor simulation 
in which pions and kaons are weighted with the corresponding probabilities of mimicking a 
muon signal. In P + HF events, the number of sequential semileptonic decay-candidates 
(29262±850) is correctly modeled by the rate of sequential decays of single 6-quarks predicted 
by the simulation (29190 ± 1236). This number is 2.5% of the P + HF total contribution 
and (6.9 ± 0.4)% of the bb contribution (424506 ± 18454 events). 

In the remaining 295481 ± 9271 ghost events, the number of additional muons in an 
angular cone with cos 6 > 0.8 around a primary muon is 49142 ± 519. In the absence of 
a simulation of ghost events, that study assumes that tracks in ghost events are a 50-50% 
mixture of pion and kaons, and estimates the number of misidentified additional muons to be 
20902 ± 284. The resulting number of unaccounted ghost events with three or more muons 
is 27970 ± 538, (9.5 ± 0.4)% of the ghost events. 

As shown by Table IIXI one half of the ghost events arise from heavy flavor production 
acquired with a misidentified muon. This type of event should contain an appreciable fraction 
of additional muons due to semileptonic decays of heavy quarks. This contribution was not 
included in the estimate of Ref. 5j. 

We estimate this contribution using Kg and K*° candidates due to it — > /i and K — > fi 
misidentification, respectively. As shown in Sec. |V] there are 32445 ± 421 and 87471 ± 2217 



Therefore, the contribution of in-flight-decays to additional muons is negligible in comparison with the 
punchthrough contribution. 
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candidates, respectively. We measure the fraction of these candidates in which at least one 
of the primary muons is accompanied by an additional muon in a cos 8 > 0.8 angular cone 
around its direction. We also estimate the contribution of fake additional muon by weighting 
all hadronic tracks that pass the additional muon selection criteria, assumed to be a 50-50% 
mixture of pions and kaons, with the corresponding misidentification probabilities jlj]. 

Figure [6] flTj) shows the invariant mass distribution of Kg (K*°) candidates when at least 
one primary muon is accompanied by an additional muon or a predicted misidentified muon. 
As previously done, we fit the K s distributions with two Gaussian functions to model the 
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FIG. 6: Invariant mass distribution of K s candidates accompanied by (left) an additional muon 
and (right) a misidentified muon. Lines represent the fits described in the text. 



signal and a straight line to model the background. The K*° distribution is fitted with a 
Breit-Wigner function plus a fourth order polynomial. 

The fits return 4572 ±91 and 1954 ± 109 events in which a K s meson is accompanied 
by an additional muon and by a fake muon, respectively. The fits return 10176 ± 739 and 
5230 ± 493 events in which a K*° meson is accompanied by an additional muon and by a 
fake muon, respectively. 

As shown in Fig. E]for events triggered by K s misidentifications, sometimes the additional 
muon is contributed by the second prong of the K s decay. Figure M shows the invariant mass 
distribution of primary and additional muons that pass the analysis selection. The usual fit 
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FIG. 7: Invariant mass distribution of K*° candidates accompanied by (left) an additional muon 
and (right) a misidentified muon. Lines represent the fits described in the text. 

yields 403 ± 33 events in which the additional muon is mimicked by the second leg of the 
Kg decay that also produced the primary muon. We remove this contribution to evaluate 
the fraction of real muons accompanying K — > \i misidentifications. We will add it for the 
fraction of events triggered by misidentified K% decays. 
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FIG. 8: Invariant mass distribution of K s candidates reconstructed using primary and additional 
muons. The line represents the fit described in the text. 
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After removing the predicted numbers of fake muons, the fraction of events with an 
identified Kg meson that contain additional real muons is (6.83 ±0.45)%. It is (5.6 ±1.05)% 
for events with an identified K*° meson. The average of the two fractions is (6.6±0.4)%. We 
multiply this fraction by the number of ghost events due to ordinary sources (241507±10284) 
to predict the number of real muons in events due to heavy flavor that are classified as ghost 
events because one of the primary muons was produced by a pion or kaon in-flight-decay. 
This procedure yields a slight overestimate because, according to the simulation tuned with 
the data, (53.4 ±0.4)% of the K — )■ \i misidentifications are due to events with heavy flavors, 
whereas (51.3 ± 0.6)% of the ghost events arise from heavy flavor production. As shown by 
Table |Xl the improved estimate still does not account for 12169 ± 1319 ghost events with 
additional real muons. 



TABLE X: Number of additional muons in ghost events are compared to the number of expected 
fake muons and real muons from heavy flavor decays. The data correspond to an integrated 
luminosity of 1426 pb _1 . 



Data 


Fakes 


K s second leg 


Heavy Flavor 


49142 ± 519 


20902 ± 284 


147 ± 15 


15924 ± 1179 



VIII. CONCLUSIONS 

This article reports an improved undestanding of the dimuon samples acquired by the 
CDF experiment. One dataset, corresponding to an integrated luminosity of 742 pb _1 , 
consists of 743006 events containing two central (|r/| < 0.7) primary (or trigger) muons, 
each with transverse momentum > 3 GeV/c, and with invariant mass larger than 5 
GeV/c 2 and smaller than 80 GeV/c 2 . These data are split into two subsets: one, referred 
to as P ± HF, consisting of 589015 ± 5074 events in which both muons originate inside 
the beam pipe of radius 1.5 cm; and one, referred to as ghost, consisting of 153991 ± 5074 
events in which at least one muon originates beyond the beam pipe. The study in Ref. J4] 
shows that the number and properties of P + HF events are correctly modeled by the 
expected contributions of semileptonic heavy flavor decays, prompt quarkonia decays, Drell- 
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Yan production, and instrumental backgrounds due to punchthrough of prompt or heavy- 
flavored hadrons which mimic a muon signal. A previous study 5] has investigated significant 
sources of ghost events, such as in-flight-decays of pions and kaons and hyperon decays. That 
study could account for approximately half of the ghost events but was unable to asses the 
uncertainty of the in-flight-decay prediction. The present study shows that the herwig 
parton-shower generator provides an accurate model of the data. The large discrepancy in 
the previous study was generated by not including the contribution of final states in which a 
6(c) hadron decays semileptonically and the second muon is produced by the in-flight-decay 
of a particle in the recoiling jet. After tuning by a few percent the pion and kaon rates 
predicted by the simulation with a fit to the data, we show that ordinary sources, mostly 
in-flight-decays, account for 125665 ± 5351 of the 153991 ± 5074 ghost events isolated in the 
sample of 743006 dimuons. 

For comparison, a DO study has used a similar dimuon sample to set a limit {3]] of 
(0.4 ± 0.26 stat ± 0.53 syst)% to the fraction of muons produced at a distance larger than 
1.6 cm from the beamline including pion and kaon in-flight-decays. This appears to be in 
contradiction with the present result, and also with a recent estimate 3jJ of the fraction of 
K — > fi and tt — > fi contributions in the DO subset of same-charge dimuons (~ 40%). 

The present study also improves a previous estimate [5] of the content of additional muons 
with pt > 2 GeV/c and \t]\ < 1.1 in ghost events. We find that (23±6)% of the unaccounted 
ghost events contain additional real muons. For comparison, the fraction of bb events that 
contain additional muons due to sequential semileptonic decays is (6.9 ± 0.4)%. 

Both results presented in this article have implications for measurements derived in 
dimuon datasets without properly accounting for the presence of ghost events. As an exam- 



ple, the measurement of the dimuon charge asymmetry performed by the DO experiment [31 1 
estimates the fraction of K — > n and 7r — > \l misidentifications with a similar method. After 
removing this background, the remaining muon pairs with same charge are attributed to bb 
production. The present study shows that, after removing this type of misidentified muons, 
the data set still contains an additional component that cannot be accounted for with ordi- 
nary sources. The size of this component, equally split in opposite and same sign pairs jsj, 
is (12.8 ± 3.2)% of the total number of dimuons due to bb production. 
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